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In  this  work  we  present  a  model  based  study  to  investigate  the  presence  of  anion  exchangers  in  direct 
methanol  fuel  cell  (DMFC)  systems.  It  is  well  known  that  environmental  or  fuel  impurities  lead  to 
accumulation  of  harmful  anions,  such  as  chloride,  in  the  system.  However,  due  to  DMFC  anodic  reaction, 
a  carbonate  system  is  present.  These  corbanate  anions  have  to  be  taken  into  account  for  the  anion 
exchanger  design  and  placement  as  well  as  for  the  system  operation  strategy  with  and  without  anion 
exchanger,  which  is  the  objective  of  this  study.  For  this  purpose,  the  expected  amount  of  harmful 
chloride  ions  in  a  DMFC  system  is  estimated,  and  that  of  carbonate  ions  is  calculated  with  a  model  of  the 
carbonate  system  in  a  DMFC  system.  The  predicition  of  durability  and  dimensions  of  an  anion  exchanger 
is  based  on  a  monovalent  anion  exchange  model.  The  design  of  gas  liquid  separators  in  the  DMFC  system 
has  a  major  influence  on  the  amount  of  dissolved  carbon  dioxide,  which  is  crucial  for  durability  and 
dimension  of  a  system  integrated  anion  exchanger.  Finally,  feasible  positions  of  anion  exchanger  in  a 
DMFC  system  are  elaborated  to  fulfill  the  needs  for  long  term  and  stable  DMFC  operation. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Direct  methanol  fuel  cell  (DMFC)  systems  are  attractive  energy 
sources  for  portable  appliances.  In  order  to  prevent  DMFC  degra¬ 
dation,  the  accumulation  of  ionic  species  in  the  DMFC  system  from 
various  sources  should  be  avoided,  as  they  negatively  affect  DMFC 
system  performance  [1-4].  Specifically,  Cl-  enters  the  system  as 
an  impurity  in  fuel  cell  components,  such  as  catalyst  layer  or 
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membrane,  or  as  a  pollutant  in  the  charged  air.  Due  to  water  re¬ 
covery  from  the  cathode  side  and  the  presence  of  NH4CI  in  the 
supplied  air  [3],  an  enrichment  of  Cl-  in  the  anode  loop  is  possible. 
Besides,  an  enrichment  of  metal  and  non-metal  ions  in  the  anode 
solution  and  in  the  membrane  electrode  assembly  (MEA)  of  the 
DMFC  was  found  to  cause  significant  performance  issues  [5]. 

To  avoid  the  before  mentioned  accumulation  of  harmful  ionic 
impurities  in  the  anode  loop,  including  anion  exchangers  into  the 
DMFC  system  has  been  proposed  [5].  However,  the  application  of 
anion  exchangers  in  DMFC  systems  is  strongly  influenced  by  the 
DMFC  system  itself.  Due  the  production  of  CO2  in  the  DMFC  anode 
reaction,  CO2  dissolves  and  forms  carbonate  species  such  as 
hydrogen  carbonate  (or  bicarbonate),  HCO3 ,  and  carbonate  ions, 
CC>3_  [6,7].  In  a  DMFC  system  with  an  anion  exchanger  in  the  anode 
loop,  these  anionic  carbon  species  rival  with  anions  from  impurities 
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for  the  ion  exchange  and  can  therefore  determine  the  durability 
and  dimension  of  the  applied  anion  exchanger. 

In  this  paper,  we  evaluate  the  feasibility  of  an  anion  exchanger 
and  its  placement  in  a  DMFC  system.  First  we  estimate  the  amount 
of  harmful  ions  like  Cl-  as  well  as  the  amount  of  carbonate  species 
which  occur  at  different  positions  in  the  DMFC  system.  This  is  done 
by  means  of  a  systematic  investigation  of  the  present  carbonate 
system  in  the  DMFC  system.  A  model  based  approach  is  applied  to 
account  for  the  amount  of  dissolved  carbonate  species  as  a  function 
of  operating  parameters  of  the  fuel  cell,  such  as  temperature  and 
methanol  concentration.  Secondly,  a  monovalent  anion  exchanger 
model  is  set  up  and  charged  with  different  flow  rates  of  solution 
containing  varying  concentrations  of  before  mentioned  anions  to 
determine  its  durability  and  dimensions.  As  a  result,  we  show  that, 
depending  on  the  gas  phase  above  the  solution,  carbonate  systems 
occur  in  the  DMFC  system:  the  semi-open  carbonate  system  is 
expected  in  the  gas  liquid  separator,  whereas  the  closed  carbonate 
system  exists  in  the  DMFC  anode  compartment  and  all  other  parts 
of  the  anode  loop,  especially  in  the  anion  exchanger.  The  occur¬ 
rence  of  different  carbonate  systems  influences  the  position  of  the 
anion  exchanger.  Based  on  these  results,  we  demonstrate  which 
anion  exchanger  positions  in  the  DMFC  system  are  feasible  to  fulfill 
the  needs  for  long  term  and  stable  operation. 

2.  Methods  and  preliminary  considerations 

In  the  following  section,  the  expected  amount  of  ion  species  of 
chloride  and  carbonate  ions  is  estimated  and  the  essential  model 
equations  for  this  work  are  derived.  The  DMFC  system  considered 
in  this  work  is  sketched  in  Fig.  1.  It  consists  of  the  DMFC,  the  gas 
liquid  separator  at  the  anode  side,  the  liquid  gas  separator  at  the 
cathode  side,  the  methanol  supply  with  methanol  tank  and  fuel 
pump,  and  the  respectively  placed  anion  exchanger  in  the  anode 
loop. 

Expected  amount  of  chloride  impurities.  Chloride  ion  impurities 
can  enter  the  DMFC  anodic  loop  via  the  neat  methanol  solution 
used,  which  possesses  a  given  chloride  impurity  of  0.5  ppm  that 
results  in  1  x  10-5  mol  l-1  Cl-  for  1  M  methanol  solution.  Another 
source  of  chloride  ions  is  the  liquid  gas  separator  behind  the  DMFC 
cathode  compartment  where  product  water  is  recycled  and  fed 
back  to  the  anodic  loop.  Due  to  the  usage  of  atmospheric  oxygen  the 
concentration  of  methyl  chloride  (CH3CI),  which  is  the  most 
abundant  halocarbon  in  the  atmosphere  and  therefore  the  most 
important  source  for  chloride  anions,  reaches  average  values  of 
approximately  500  ppt  [3  .  Due  to  the  very  low  solubility  of  chlo¬ 
ride  in  water  [8  ,  we  estimate  a  concentration  of  2  x  10“7  mol  l-1 
dissolved  Cl-  for  this  work. 


Fig.  1.  Schematic  of  the  DMFC  system  with  different  options  for  positioning  the  anion 
exchanger.  The  DMFC  system  is  described  in  more  detail  elsewhere  [23]. 


2.1.  Expected  amount  of  carbonate  ion  species 

If  gaseous  CO2  is  in  contact  with  an  aqueous  solution,  it  partly 
dissolves  in  the  solution  and  a  carbonate  system  evolves.  In  the 
following,  we  give  a  brief  introduction  on  the  carbonate  system  and 
investigate  its  appearance  in  DMFC  systems. 

Dissolved  CO2  yields  three  different  inorganic  species  in 
aqueous  solutions:  dissolved  carbon  dioxide,  C02(diss);  bicarbon¬ 
ate,  FICO3 ;  and  carbonate  ion,  CO3- .  All  together  constitute  the 
carbonate  system  [6]  and  their  individual  concentrations  are 


related  as  follows: 

[C02(diss)]  =H£02  -p'Co2  (1) 

[h+]  -  [HCOj]  =  [C02(diss)]  Xal  (2) 

[h+]  •  [cof  ]  =  [HCO 3]-I<a2  (3) 

[H+]  •  [OH"]  =  Kw  (4) 


where  concentrations  in  liquid  phase  are  indicated  with  square 
brackets  and  the  unit  of  H£02 ,  Flenry’s  law  coefficient,  is  mol  l-1  [9]. 
I<  values  indicate  the  pressure  and  temperature  dependent  disso¬ 
ciation  constants.  JCW  is  referred  to  as  the  ion  product  of  water.  The 
respective  constants  used  in  the  present  work  are  given  in  Table  1. 

Consequently,  the  carbonate  system  is  described  completely 
with  the  use  of  Equations  (l)-(4),  the  sum  of  all  dissolved  carbon 
dioxide  species  (DIC)  (5)  and  the  electroneutrality  condition  (6): 

DIC  =  [HCOj]  +  [CO^-]  +  [C02(diss)]  (5) 

0  =  [H+]  ~~  [hc03 ]  “2  -  [COI-]  -  [0H“]  (6) 

The  influence  of  methanol  on  the  dissolution  of  carbon  dioxide 
is  given  by  Ref.  [13].  Since,  a  highly  diluted  methanol  water  solution 
is  expected  at  the  anode  outlet,  the  influence  of  methanol  on  the 
carbonate  system  can  be  neglected  for  DMFC  operation.  This 
accompanied  by  considering  the  methanol  influence  on  Flenry’s 
law  coefficient  for  the  carbonate  system  as  described  in  Section  2.2. 

Depending  on  the  properties  or  existence  of  a  gas  phase  above  a 
solution,  the  carbonate  system  can  be  classified  as  open  (a),  closed 
(b)  or  semi-open  (c)  carbonate  system  [6].  The  open  system  is  the 
most  complex,  as  found  e.g.  in  sea  water  systems  in  contact  with 
environmental  air.  We  consider  it  as  not  relevant  for  DMFC  oper¬ 
ation.  The  semi-open  system  possesses  a  gas  phase  with  a  constant 
amount  of  gaseous  carbon  dioxide  with  a  constant  partial  pressure 
p'C02  as  it  may  occur  in  certain  parts  of  the  DMFC  system.  The  closed 
system  possesses  no  gas  phase  above  the  solution  that  interacts 
with  the  environmental  atmosphere.  A  schematic  representation  of 
the  carbonate  systems  used  in  this  work,  is  depicted  in  Fig.  2.  In  the 
same  figure,  an  overview  on  the  expected  amounts  of  carbonate 
species  is  given  for  the  closed  carbonate  system  and  the  semi-open 
carbonate  system.  The  minimal  depicted  partial  pressure  is  set  to 
the  expected  partial  pressure  of  carbon  dioxide  in  the  atmosphere 
of  about  3.535  x  10  5  MPa,  whereas  the  maximal  depicted  partial 
pressure  is  equal  to  atmospheric  pressure  of  0.101325  MPa  at  298  K. 
According  to  p'C02  =  0.101325  MPa  given  by  the  anode  side  gas 
liquid  separator,  the  aqueous  methanol  solution  will  be  present 
with  a  maximal  possible  concentration  of  dissolved  carbonate  of 
cco2(diss)  =  2.26  x  10_2mol  l-1.  Table  2  lists  all  relevant  DMFC 
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Table  1 

Equilibrium  constants  used  in  this  work  reference. 


Constant 

Expression 

Unit 

Reference 

Hco2 

exp(  6.8346  12f17  3766800  +  299700000^) 

MPa 

[10] 

/Cal 

exp(  356.3094  0.06091964  T  +  21834- 37  +  126.8339  log(T)  1684915) 

mol  l-1 

[11] 

/Ca2 

exp(  107.8871  0.032528493  T  +  515f 79  +  38.92561  log(T)  563£3-9) 

mol  l-1 

[11] 

/<w 

exp (2 1.6470  42394000  3.7546  log(T)) 

mol2  r2 

[11],  respectively  [12] 

components,  indicated  with  the  respective  carbonate  systems  as 
they  might  appear  in  DMFC  systems. 

2.2.  Model  of  DMFC  anode  compartment 


Knowing  that  methanol  is  fed  with  excess  by  the  factor  of  Amc 
the  molar  flow  rate  of  methanol  entering  the  anode  is: 

^Me,in  —  ^Me'^Me,need  (8) 


In  the  following,  a  steady  state  model  for  the  DMFC  anode  loop,  as 
illustrated  in  Fig.  1,  containing  an  anion  exchanger  at  different  posi¬ 
tions  is  presented.  It  consists  of  a  DMFC  anode  compartment  model 
and  models  for  the  other  anodic  loops  parts  to  describe  the  expected 
amount  of  carbon  species  in  each  part  depending  on  DMFC  operating 
conditions.  The  DMFC  anode  compartment  model  consists  of  mass 
balances  describing  the  liquid  and  gaseous  species. 

It  is  assumed  that  the  considered  DMFC  operates  at  a  temper¬ 
ature  Tfc  and  atmospheric  pressure  patm.  The  molar  flow  rate  of 
produced  carbon  dioxide  is  equal  to  the  molar  flow  rate  of  meth¬ 
anol  needed  for  the  DMFC  anode  reaction.  According  to  Faraday’s 
law,  it  follows: 


^Me,need  — 


^  *  ^ceM 

6F 


^C02,prod 


(7) 


Therefore,  the  volumetric  flow  rate  of  fed  methanol  from  the 
fuel  pump  is  expressed  by: 


^Me  = 


^Me,need '  ^Me 
PMe(^cool) 


(9) 


The  density  of  methanol  pMe  is  determined  with  data  by  Lide 

[14]. 

The  molar  flow  rate  of  solution  hfeed  entering  the  anode  is  the 
sum  of  the  molar  flow  rates  of  methanol  and  water: 


^feed  —  ^Me,in  +  ^H20,in  (10) 

Due  to  the  fact  that  the  system  is  fed  with  a  diluted  water 
methanol  solution  (<3  wt  %),  the  volumetric  flow  rate  of  solution 
entering  the  anode  Ffeed  is  approximated  by: 


10“ 4  10” 3  10” 2  10” 1  300  310  320  330  340  350 


Pco2  lMPal 


T  [K] 


Fig.  2.  The  carbonate  system  and  a)  its  classification  according  to  Stumm  and  Morgan  [6];  b)  the  concentration  of  the  ionic  species  for  a  closed  carbonate  system  at  1  atm  and  298  K 
and  for  semi-open  conditions:  c)  versus  partial  pressure  range  and  d)  versus  temperature  range. 
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Table  2 

Assignment  of  classification  of  the  carbonate  system  to  relevant  DMFC 
components. 


Component 

Type 

DMFC  (anode  compartment) 

Closed 

Cooler 

Closed 

Separator 

Semi  open 

Mixer 

Closed 

Circulation  Pump 

Closed 

_  Hfeed  *  ^feed  ^  ^feed  '  ^H20  Mil 

feed  "  PfeedTcool)  PH2oTcool)  1  ; 

where  it  is  assumed  that  the  density  and  the  molar  mass  of  the 
diluted  methanol  water  solution  possesses  about  the  same  density 
as  pure  water  [15].  The  density  of  water  is  calculated  with  an 
empirical  equation  [16]. 

Furthermore,  the  molar  amount  of  solution  lost  in  the  anode 
compartment  is  the  sum  of  lost  methanol  and  water  due  to  cross¬ 
over  from  anode  to  cathode  compartment  and  anode  reaction 
(subscript  xover  and  reaction,  respectively),  and  the  molar  flow  rate 
of  solution  exiting  the  anode  compartment  is  denoted  as  hout. 
Considering  this,  it  follows: 

ftout  —  ftfeed  —  ^Me, reaction  —  ^Me,  xover  —  ^H20, reaction  —  ^H20,  xover 

(12) 

Thereby,  the  terms  of  reaction  for  water  and  methanol  are 
stoichiometrically  expressed  by  Equation  (8).  Water  crossover  is 
accounted  for  with  a  water  drag  coefficient  k  [17]: 

^H20, xover  =  k-  (13) 


described  separately,  but  determines  the  operating  temperature  for 
the  gas  liquid  separator  with  its  temperature  rcooi. 

The  dissociation  reactions  in  the  carbonate  system  (Equations 
(2)— (4))  take  place  in  the  range  of  10-2  s  and  10  s  [18].  The  slowest 
process  is  the  dissolution,  Equation  (1)  [6].  We  assume  the  resi¬ 
dence  time  of  the  methanol  water  solution  in  all  system  compo¬ 
nents  of  the  given  standard  DMFC  system  to  be  sufficient  to  obtain 
thermodynamic  equilibrium  regarding  the  carbonate  system. 

2.3.  Possible  positions  of  the  anion  exchanger 

The  DMFC  system  considered  in  this  work,  as  sketched  in  Fig.  1, 
operates  at  constant  current  and  in  steady  state.  This  implies  that  a 
constant  volumetric  flow  of  solution  enters  and  leaves  each 
component  in  the  anode  loop,  e.g.  the  gas  liquid  separator.  Further¬ 
more,  the  gas  liquid  solution  is  assumed  to  be  in  thermodynamic 
equilibrium  regarding  the  carbon  dioxide  species  with  its  dissolved 
species. 

At  position  1,  the  anion  exchanger  is  investigated  at  the  meth¬ 
anol  tank  outlet  fed  with  constant  chloride  impurity  and  a  volu¬ 
metric  flow  rate  FMe  according  to  a  set  DMFC  system  steady  state. 
The  given  chloride  impurity  by  the  manufacturer  of  0.5  ppm  results 
in  about  1  x  10-5  mol  l-1  Cl-  in  the  pure  methanol  solution.  For  the 
anion  exchanger  at  position  2,  the  water  regained  from  the  cathode 
side  is  saturated  with  CO2  with  environmental  partial  pressure  of 
carbon  dioxide  of  3.535  x  10-5  MPa.  This  results  in  a  concentration 
of  [C02(diss)]  ~  1  x  10-5  mol  l-1  in  the  recycled  water.  At  position 
3,  the  anion  exchanger  is  considered  to  be  at  the  anode  compart¬ 
ment  inlet  fed  with  Ffeed  containing  dissolved  carbon  dioxide 
cco2(diss)  Cl-.  An  anion  exchanger  at  this  position  has  the  task  to 
remove  all  harmful  anionic  impurities  completely.  However,  the 
expected  concentration  of  carbonate  species  anions  is  marginally 
higher  then  for  example  the  Cl-  concentration  at  this  position. 


Methanol  influences  the  equilibrium  constant  of  the  gas  liquid 
interaction  between  carbon  dioxide  and  the  solution,  and  it  impacts 
the  dissociation  of  the  carbonate  system  species.  A  vast  set  of 
experimental  data  for  Henry’s  law  constants  in  pure  water  systems 
was  previously  investigated  [10,11],  presenting  temperature 
dependent  correlations  valid  for  atmospheric  pressure.  The  impact 
of  methanol  concentration  on  Henry’s  law  constants  for  methanol 
water  solutions  with  CO2  is  given  in  literature  [13].  A  comparison  of 
Henry’s  law  constants  with  and  without  methanol  influence  is 
given  in  Appendix  B.  Henry’s  law  constants  from  the  above 
mentioned  authors  are  plotted  as  function  of  temperature  in 
Fig.  B.l.  It  has  to  be  noted  that  the  given  Henry’s  law  constants  are 
converted  to  MPa  units  to  ensure  comparability  [9]. 

Methanol  crossover  is  neglected  for  this  study,  as  the  difference 
in  concentration  will  not  effect  the  solubility  of  CO2,  compare 
Fig.  B.l,  so  that  hMe  xover  =  0.  Finally,  this  results  in: 

nout*fe-2-kV^f“  (14) 

\*Me  J  b  •  r 

where  XMe  represents  the  molar  fraction  of  methanol  in  the  liquid 
phase  which  depends  on  the  molar  concentration  of  methanol  and 
water  in  the  anode  feed  solution. 

Assuming  the  same  conditions  as  in  Equation  (11),  the  volu¬ 
metric  flow  rate  of  the  liquid  phase  is  approximated  by: 


Tout  ~ 


n0uf^H2o 

PH2o(Tcool) 


(15) 


2.4.  Anion  exchanger  model 

The  considered  anion  exchanger  is  a  column  filled  with  anion 
exchange  resin  particles.  In  the  following  we  assume  an  ideal 
monovalent  anion  exchange  of  two  kinds  of  anions:  either  HCO3 , 
originating  from  the  carbonate  system,  or  Cl-,  exemplary  for  impu¬ 
rities.  The  anion  exchange  column  is  fed  with  a  volumetric  flow  rate 
from  the  DMFC  system  containing  the  anions  A  and  B,  where  A  is  the 
anion  which  is  intended  to  be  exchanged,  as  HCO3  or  CO2-  or  Cl-, 
respectively.  Anion  B  is  considered  to  be  OH-  and  is  already  loaded 
onto  the  resin  particles.  Furthermore,  the  packed  resin  particles  are 
equally  distributed  possessing  a  uniform  particle  size.  Mass  transfer 
resistance  in  the  macro  pores  of  the  resin  and  axial  dispersion  in  the 
packed  column  as  well  as  possible  back  flows  inside  the  column  are 
neglected. 

In  a  DMFC,  the  present  ionic  species  is  HCO3  rather  than  CO2-. 
We  therefore  consider  solely  HCO3  as  the  ion  which  will  be 
intentionally  removed  from  a  solution  passing  the  anion  exchanger. 

If  the  anion  exchange  resin  particles  are  completely  loaded  with 
ion  HCO3  or  CO2-  or  Cl-,  respectively,  it  is  considered  to  be  full  and 
the  durability  of  the  anion  exchanger  is  reached. 

To  determine  the  durability  of  the  anion  exchanger,  a  model  for 
anion  exchange  (AIXC)  columns  is  formulated.  It  combines  an  ideal 
plug  flow  model  with  the  estimation  of  a  local  anion  exchange 
equilibrium  according  to  an  equilibrium  isotherm  [19].  The  anion 
exchanger  durability  is  the  ratio  of  column  dimensions  and  prop¬ 
agation  velocity  of  the  ion  species  A  inside  the  column: 


With  Equations  (7)  and  (14),  the  molar  flow  of  gas  and  liquid  L-J-D2  (*+(1-*).^)  ,  x 

exiting  the  anode  and  entering  the  cooler,  and  consequently  the  gas  £aixc  = - - — p - —  (16) 

liquid  separator,  can  be  described.  The  anode  cooler  is  not  IE 
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Table  3 

Reference  operating  conditions  at  steady  state  for  DMFC  operation  and  anion 
exchanger  parameter. 


Operating  parameter 

Value 

Unit 

Patm  =  Pc  o2 

0.101325 

MPa 

Tpc 

333 

K 

Trool 

310 

K 

Ncell 

1 

- 

/ 

3.5 

A 

^Me 

3 

- 

Qvie 

1 

mol  1 

k 

0 

- 

D 

4 

mm 

L 

40 

mm 

e 

0.4 

- 

Q 

0.8 

mol  1 

The  inter  particle  porosity  of  the  anion  exchanger  resin  in  Equation 
(16)  is  expressed  by  e  and  constant  volumetric  flow  and  concentration 
CA,in  and  a  low  concentration  cb, in  in  the  fed  solution  are  assumed. 

All  subsequently  presented  calculations  for  the  anion  exchanger 
durability  are  calculated  with  the  chosen  reference  operating 
conditions  given  in  Table  3. 

3.  Results  and  discussion 

In  the  following  section,  we  present  the  results  of  the  model 
based  analysis  and  elucidate  the  anion  exchanger  lifetime  accord¬ 
ing  to  its  placement  in  the  DMFC  system.  It  was  shown  in  the 
previous  section  that  the  concentration  of  carbonate  species  is 
much  higher  then  the  concentration  of  impurities  in  the  system. 
Therefore,  we  first  estimate  the  durability  of  the  anion  exchanger 
solely  for  cr  and  then  solely  for  the  presence  of  HCO3 . 

3.1.  Anion  exchanger  durability  based  on  chloride  impurities 

For  the  three  different  possible  positions  of  the  anion  exchanger  in 
the  DMFC  anode  loop  the  resulting  anion  exchanger  lifetimes  are 
given  in  Table  4.  At  position  1,  the  anion  exchanger  is  investigated  at 
the  methanol  tank  outlet  fed  with  FMe  =  15.04  pi  min-1  at  chosen 
reference  conditions  and  a  low  chloride  concentration,  this  results  in 
high  anion  exchanger  durability.  At  position  2,  behind  the  liquid  gas 
separator  on  the  cathode  side,  we  assume  the  dissolved  amount  of 
chloride  to  be  much  lower  than  the  amount  of  dissolved  carbon  di¬ 
oxide  due  to  the  very  low  solubility  of  chloride  in  water  [8].  Conse¬ 
quently,  while  Cl-  will  also  be  removed,  solely  carbonic  anions  are 
considered  to  determine  the  lifetime  of  an  anion  exchanger  at  this 
position  (see  next  subsection).  At  position  3,  the  anion  exchanger  is 
considered  to  be  at  the  anode  compartment  inlet  fed  with  Ffeed  con¬ 
taining  dissolved  carbon  dioxide  cC02(diss)  and  Cl-.  An  anion  exchanger 
at  this  position  has  the  task  to  remove  all  anionic  impurities 
completely.  Assuming  that  the  anions  are  caused  by  the  impurities  of 
the  refreshed  methanol  solution,  the  Cl-  concentration  is  approxi¬ 
mated  as  2  x  10  7  mol  l-1.  This  results  in  an  anion  exchanger  durability 
for  chloride  of  about  3  years,  which  is  a  suitable  lifetime  compared  to 
other  DMFC  system  components. 


Table  4 

Calculated  anion  exchanger  lifetimes  depending  on  the  concentration  of  chloride 
ions  at  the  corresponding  position  in  the  DMFC  anode  loop. 


Position 

Expected  concentration  mol  1  1 

Anion  exchanger  lifetime 

1) 

1  x  lO”5 

2.5  years 

2) 

~0 

OO 

3) 

2  x  lO"7 

3  years 

3.2.  Anion  exchanger  durability  based  on  carbonate  ions 

In  addition  to  the  chloride  impurities  also  carbonate  ions  are  pre¬ 
sent  in  the  anode  loop  due  to  the  reaction  of  methanol  and  water  in  the 
anode  compartment  of  the  DMFC.  Depending  on  the  position  in  the 
loop,  the  anion  exchanger  will  be  fed  with  a  different  amount  of  car¬ 
bonate  ions.  For  the  positions  given  in  Fig.  1,  the  resulting  anion 
exchanger  durability  are  presented  in  Table  5.  At  position  1,  the  anion 
exchanger  is  investigated  at  the  methanol  tank  outlet,  where  we 
expect  only  pure  methanol  without  carbonates  to  be  applied.  Conse¬ 
quently  the  carbonate  ion  concentration  and  other  impurities  than 
chloride  in  the  methanol  solution  are  neglected  at  this  position.  For  the 
anion  exchanger  at  position  2,  the  water  regained  from  the  cathode 
side  is  saturated  with  CO2  with  environmental  partial  pressure  of 
carbon  dioxide  of 3.535  x  1CT5  MPa.  In  real  DMFC  operation  this  value 
can  even  be  higher  due  to  the  methanol  crossover,  as  CO2  is  produced 
in  a  side-reaction  [20]  and  the  partial  pressure  will  be  higher  than  the 
environmental  partial  pressure. 

The  volumetric  flow  rate  of  the  recycled  water  is  given  by  the 
reaction  (nH2o, reaction)  as  well as  by  the  water  crossover  (hH2o,xover)> 
and  with  k  =  10  it  is  7.30  pi  min-1.  Therefore  the  durability  of  the 
anion  exchanger  is  0.69  years,  if  all  ionic  carbon  species  are 
removed.  The  long  durability  is  mainly  due  to  the  very  low  flow  rate 
of  recycled  water  and  can  be  considered  feasible  for  portable  DMFC 
systems.  At  position  3,  the  anion  exchanger  is  considered  to  be  at 
the  anode  compartment  inlet  fed  with  Ffeed  =  1.10  ml  min-1  at 
standard  operating  conditions.  According  to  p'C02  =  0.101325  MPa 
given  by  the  anode  side  gas  liquid  separator,  the  methanol  solution 
contains  the  maximal  possible  concentration  of  cC02(diss)  =  2.26  x 
10~2 3mol  l-1 .  This  results  in  an  anion  exchanger  durability  of  about 
only  0.144  h.  Since,  dissolved  carbon  dioxide  has  the  main  influence 
on  the  anion  exchanger  durability  at  position  3,  durability  might 
not  be  suitable  for  portable  DMFC  system  operation. 

Another  way  to  influence  the  amount  of  dissolved  carbon  dioxide 
is  given  by  the  design  of  the  gas  liquid  separator  in  the  anode  loop  of 
the  DMFC  system.  The  H+  concentration  and  the  pH  value  in  the 
anode  loop  are  mainly  determined  by  the  partial  pressure  of  carbon 
dioxide  in  the  gas  phase.  To  illustrate  the  latter,  Fig.  3  shows  the  pH 
value  in  the  semi-open  carbonate  system  as  a  function  of  system 
temperature  and  partial  pressure  of  carbon  dioxide.  Thereby,  the 
minimal  pH  value  is  3.90,  depicted  in  the  bottom  corner  at  maximal 
partial  pressure  and  minimal  system  temperature.  By  implication, 
pH  =  5.75  in  the  upper  corner  is  assigned  to  minimal  partial  pressure 
and  maximal  system  temperature.  The  pH  value  at  constant  partial 
pressure  changes  only  slightly  as  function  of  temperature,  for 
example  atp^  =  0.101325  MPa  from  pH  3.9060  to  4.0940. 

This  implies  for  the  cooler  and  the  gas  liquid  separator  operation 
that  a  change  in  operating  temperature  does  not  change  the  pH 
value  significantly.  Consequently,  the  design  of  a  gas  liquid  sepa¬ 
rator  in  a  DMFC  (compare  Fig.  3  1)  and  2))  might  be  considered 
solely  regarding  the  pressure  of  carbon  dioxide.  Assuming  an 
application  of  a  micro  separation  unit,  as  presented  by  Kraus  et  al. 
[21],  a  significant  higher  pH  value,  which  corresponds  to  a  lower 
concentration  of  dissolved  carbonate  ions,  is  expected  at  the  unit 
outlet.  This  is  due  to  the  fact  that  the  liquid  phase  in  this  unit  has 
direct  contact  with  the  environmental  air.  Therefore  the  design  of 

Table  5 

Calculated  anion  exchanger  durability  depending  on  the  concentration  of  carbonate 
ions  at  the  corresponding  position  in  the  DMFC  anode  loop. 

Position  Expected  concentration  mol  1_1  Anion  exchanger  durability 

1)  =0 

2)  =1  x  10'5  0.69  years 

3)  2.26  x  lO"2  0.144  h 
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Fig.  3.  pH  value  in  the  semi-open  carbonate  system  as  function  of  temperature  and  partial  pressure  of  carbon  dioxide.  The  design  of  the  semi-open,  capillary-based  gas  liquid 
separators  [21]  with  1)  an  open  separation  area,  where  the  degassed  methanol  solution  stays  in  contact  with  environmental  air,  and  2)  a  covered  separation  area,  where  C02  gas  is 
accumulated  up  to  the  environmental  pressure  of  1  bar. 


micro  channel  based  separation  units  can  have  a  major  influence  on 
the  amount  of  the  carbonic  anion  species  in  DMFC  systems. 

3.3.  Dimension  of  anion  exchanger  depending  on  system  and  design 

According  to  Equation  (16),  clear  trends  to  specify  the  anion 
exchanger  durability  are  observed.  Increasing  the  column  length 
L  and  the  anion  exchange  capacity  Q.  increases  the  anion  exchanger 
durability  linearly.  Extending  the  column  diameter  D  increases  the 
anion  exchanger  durability  quadratically.  An  increase  in  inter 


particle  porosity  e  decreases  the  anion  exchanger  durability  linearly. 
Since  portable  operation  favors  a  reduction  of  column  dimensions, 
the  column  diameter  D  is  an  appropriate  way  to  scale  the  anion 
exchanger  durability.  Fig.  4a)-d)  shows  the  influence  of  selected 
DMFC  system  and  anion  exchanger  design  parameters  on  the  pre¬ 
dicted  anion  exchanger  durability,  when  it  is  placed  directly  in  front 
of  the  DMFC  anode  inlet  at  position  2.  Thereby  all  parameters  are 
set  to  the  chosen  reference  conditions  from  Table  3,  except  the 
respective  investigated  parameter.  Fig.  4a)  to  b)  depict  the  results  for 
parameters  which  are  related  to  the  flow  rate  in  the  anode  loop.  I  and 


0  5  10  15  20  25 

k 


D  [mm] 


Fig.  4.  Anion  exchanger  durability  tret  as  function  of:  a)  fuel  cell  current,  b)  methanol  excess  ratio  AMe.  c)  water  drag  coefficient  k  and  d)  column  diameter  D. 
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X  both  increase  with  Ffeed  and  as  such  show  reciprocal  influence  on 
the  anion  exchanger  durability. 

As  presented  in  Fig.  4c),  the  water  drag  k  invokes  only  a  slight 
and  almost  linear  influence  on  anion  exchanger  durability.  This  can 
be  explained  by  water  moving  to  the  cathode  which  decreases  hou t, 
and  as  such  the  total  amount  of  carbon  dioxide  which  is  taken  up  by 
the  solution.  This  solution  is  then  diluted  in  the  mixer  with 
recovered  water  from  the  cathode  condenser  which  contains  a 
significant  lower  amount  of  ionic  carbon  species. 

Fig.  4d)  depicts  exemplary  one  result  for  anion  exchange  column 
design  parameters,  the  column  diameter.  Compared  to  the  impact 
of  DMFC  operating  system  parameters,  the  column  diameter  in¬ 
fluences  the  anion  exchanger  durability  stronger. 

4.  Conclusions 

In  this  paper,  model  based  investigations  of  the  carbonate  sys¬ 
tem  and  anion  exchangers  in  a  DMFC  system  are  presented.  It  is 
shown,  that  anions  are  not  exclusively  present  in  the  DMFC  system 
due  to  environmental  or  methanol  impurities  but  rather  by  the 
electrochemical  CO2  production  at  the  anode  compartment  and  the 
consequently  evolving  carbonate  system.  For  the  removal  of  anions, 
the  usage  of  an  anion  exchanger  at  different  positions  and  design  in 
the  DMFC  system  is  considered. 

If  the  accumulation  of  harmful  anions  like  chloride  in  the  anode 
loop  and  consequently  DMFC  degradation  is  occurring,  anion  ex¬ 
changers  are  a  suitable  tool  to  implement  in  DMFC  systems.  We  have 
shown  by  a  model  based  analysis  that  only  certain  positions  in  the 
DMFC  system  are  practical  for  example  for  portable  operation:  a  po¬ 
sition  directly  in  front  of  the  anode  inlet  is  not  suitable  for  portable 
applications  because  carbonate  ions  outnumber  any  impurities  by  far 
and  as  such  lead  to  very  low  anion  exchanger  durability.  The  position 
behind  the  methanol  tank  as  well  as  the  position  in  the  recycling  line 
offer  suitable  lifetimes  and  anion  exchanger  sizes  for  a  portable  usage. 

Regarding  the  carbonate  system  in  a  DMFC,  the  methanol  in¬ 
fluence  can  be  neglected  within  the  expected  operation  range. 
Furthermore,  the  concentrations  of  carbonate  species  and  pFI 
values  are  mainly  influenced  by  the  partial  pressure  of  carbon  di¬ 
oxide  the  solution  has  been  exposed  to,  rather  than  by  the  system 
temperature.  Hence,  gas  liquid  separators  are  key  factors  in  DMFC 
systems  regarding  the  carbonate  ion  concentration  and  the  water 
balance  of  portable  DMFC  systems.  This  complements  the  above 
mentioned  findings  [21,22  ,  where  separators  are  considered  as  key 
parts  of  DMFC  system  design.  To  achieve  long  term  and  stable 
operation,  considering  carefully  the  gas  liquid  separators  as  well  as 
anion  exchangers  will  be  necessary  for  DMFC  system  design. 

Appendix  A.  Calculation  of  the  semi-open  and  closed 
carbonate  system 

Appendix  AA  Semi-open  system 

To  determine  the  concentration  of  dissolved  carbon  dioxide 
[C02(diss)]  in  the  semi-open  system,  Henry’s  law  is  combined  with 
the  respective  law  of  mass  action.  Thereby,  Henry’s  law  is  valid  for  a 
system  pressure  below  0.5  MPa  and  assumes  ideal  dilute  solution 
behavior.  Usually,  Henry’s  law  constant  H  is  defined  as  the  ratio  of 
the  concentration  of  a  species  in  the  solution  phase  c  and  the  partial 
pressure  p'  of  the  species  in  the  gas  phase: 


With  this  definition  the  unit  of  H  is  mol  1  1  MPa  l  However,  values  for 
H  are  expressed  in  literature  in  several  ways  using  different  units 


according  to  modified  equations  of  Henry’s  law.  Sander  [9]  gives  a 
detailed  review  on  used  Henry’s  law  constants  and  their  unit 
conversions. 

To  calculate  the  species  concentrations  the  law  of  mass  action  is 
applied  to  the  equilibrium  reactions,  which  leads  to  Equations  (1)— 
(4).  With  Equations  (5)  and  (6)  in  mind,  six  equations  for  the  seven 
unknown  concentrations  [C02(g)],  [C02(diss)],  [HCO3],  [CO§— ],  [H+], 
[OH-]  and  DIC  with  given  temperature  and  system  pressure  are 
obtained.  As  a  consequence,  one  unknown  is  additionally  needed  to 
calculate  the  remaining  species  concentrations.  Usually,  [C02(g)]  is 
given  by  the  partial  pressure  of  carbon  dioxide  p'C02  in  the  gas  phase 
above  the  solution  of  about  3.535  x  10-5  MPa  in  the  atmosphere 
according  to  Ref.  [7  .  Rearranging  the  set  of  equations,  yields  to  a 
third  order  equation  for  [H+]  as  a  function  of  the  known  value  p'C02 
and  the  equilibrium  constants: 

0  =  [h+]  3  -  (Kw  +  Kai  Hco2  -Pc o2)  ■  [h+] 

-2  Kal  ■I<a2  HCo2  p'co2  (A.2) 

whereas  it  is  important  that  Henry’s  law  constant  is  denoted  with 
mol  l-1  MPa-1  and  p'C02  with  MPa.  Only  the  one  non-complex  so¬ 
lution  of  the  three  solutions  of  Equation  (A.2)  is  physically  mean¬ 
ingful  and  yields  the  pH  value  as  the  concentration  of  H+. 
Consequently,  the  semi-open  system  can  be  validated  by  measuring 
the  pH  value  at  given  temperature  and  system  pressure. 

Appendix  A.2  Closed  system 

No  gas  phase  is  present  and  the  reaction  (A.3)  is  not  existent  in 
this  system. 

Hco2 

C02(g)^C02(diss)  (A.3) 

Consequently,  five  equations  for  the  six  unknown  concentra- 
tions  [C02(diss)],  [HCOJ],  [COf-],  [H+],  [OH  ]  and  DIC  with  given 
temperature  and  system  pressure  are  remaining.  Again,  one  un¬ 
known  is  additionally  needed  to  calculate  the  remaining  species 
concentrations.  Thereby,  the  unknown  concentration  can  be  seen 
as  entering  the  system  inlet  determining  the  other  species  con¬ 
centrations  according  to  the  thermodynamic  equilibrium. 

There  are  six  possibilities  to  calculate  the  carbon  system  species 
concentrations  in  a  closed  system  [7].  As  example,  the  result  for 
calculating  with  known  DIC  is  shown.  Combining  the  equations 
mentioned  above  leads  to  [H+]  as  a  function  of  DIC  and  the  equi¬ 
librium  constants: 

0  =  [h+]  4  +  Kal  •  [h+]  3  +  (Kal  I<a2  -  Kw  -  Kai  •  DIC)  •  [h+]  2 
+  (  -  2  Kal  Ka2  ■  DIC  -  Kal  -Kw)  •  [h+]  -  Kal  Ka2  Kw 

(A.4) 

In  general,  fourth  order  equations  like  Equation  (A.4)  have  four 
solutions.  Only  one  of  it  is  physically  meaningful.  With  it,  the 
remaining  species  concentrations  can  be  determined  according  to  the 
set  of  equations.  Again,  the  concentration  of  H+  yields  the  pH  value 
and  consequently  the  closed  system  can  be  determined  by  measuring 
the  pH  value.  In  the  limiting  case  of  DIC  becoming  zero,  the  pH  is 
approximately  7,  hence  only  water  is  present  in  the  system.  Besides, 
this  case  is  equal  to  a  semi-open  system  with  p'C02  becoming  zero. 

Appendix  B.  Methanol  in  the  carbonate  system 

We  elucidate  the  influence  of  methanol  on  the  evolved  car¬ 
bonate  system  in  a  DMFC  and  consider  the  amount  of  carbonate 
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species  which  can  be  expected  during  DMFC  operation.  Methanol 
influences  the  equilibrium  constant  of  the  gas  liquid  interaction 
between  carbon  dioxide  and  the  solution,  and  it  impacts  the 
dissociation  of  the  carbonate  system  species.  Therefore,  a  com¬ 
parison  of  Henry’s  law  constants  with  and  without  methanol  in¬ 
fluence  will  be  given  subsequently. 

Refs.  [11,10  investigated  a  vast  set  of  experimental  data  for  Hen¬ 
ry’s  law  constants  in  pure  water  systems.  They  present  temperature 
dependent  correlations  valid  for  atmospheric  pressure.  Ref.  [13] 
experimentally  determined  the  impact  of  methanol  concentration 
on  Henry’s  law  constants  for  methanol  water  solutions  with  CO2. 

According  to  Fig.  B.l  the  obtained  curves  from  Refs.  [11,10],  Henry’s 
law  constant  decreases  linearly  with  decreasing  temperatures  below 
the  maximal  expected  temperatures  in  the  DMFC  system  of  353  K. 
Besides,  it  can  be  seen  that  methanol  decreases  Henry’s  law  constant 
which  consequently  enhances  CO2  dissolution  according  to  Equation 
(1).  Methanol  molar  fractions  of  not  more  than  XMe  ~  0.02  and  tem¬ 
peratures  below  353  I<  are  expected  for  the  here  considered  portable 
autonomous  DMFC  system.  No  significant  difference  of  Henry’s  law 
constants  is  observed  in  this  region.  This  implies  that  methanol  in¬ 
fluence  on  the  carbonate  system  in  a  DMFC  system  can  be  neglected 
for  simulation  and  the  constants  for  pure  water  are  used  in  this  work. 


Fig.  B.l.  Comparison  of  Henry’s  law  constant  HCq2  as  function  of  temperature  from 
T  =  310  K  to  400  K  by  Ref.  [11  ]  (light  gray  line)  and  [10]  (dark  gray  line),  together  with 
the  influence  of  methanol  molar  fraction  with  data  by  Ref.  [13]  (black  square,  triangle 
and  circle  marker,  respectively). 


Nomenclature 

Latin  symbols 

c  concentration  (mol  l-1) 

D  column  diameter  (m) 

DIC  total  dissolved  inorganic  carbon  concentration  (mol  l-1) 
F  volume  flow  (m3  s-1) 

F  Faraday  constant  (96485.309(23)  A  s  mol-1) 

H  Henry’s  law  constant,  unit  depending  on  used  definition 

I  current  (A) 

k  water  drag  coefficient  (-) 

Kai  equilibrium  constant  (mol  l-1) 

/<a2  equilibrium  constant  (moll-1) 

/<w  equilibrium  constant  for  ion  product  of  water  (mol2  l-2) 
L  column  length  (m) 

Nceii  number  of  fuel  cells  in  a  stack  (-) 


h  molar  flow  (mol  s1) 

M  molar  mass  kg  mol-1 

p  pressure  (MPa) 

p'  partial  pressure  (MPa) 

q  anion  concentration  in  the  resin  phase  (mol  l-1) 

Q.  anion  exchange  capacity  (moll-1) 

t  time  (s) 

T  temperature  (I<) 

Greek  symbols 

e  inter  particle  porosity  (-) 

X  stoichiometry  factor  (-) 

p  density  (kg  m-3) 

Subscripts 

AIXC  anion  exchanger 

atm  physical  value  at  atmospheric  conditions 

al  regarding  the  dissociation  of  dissolved  carbon  dioxide 

a2  regarding  the  dissociation  of  hydrogen  carbonate 

cell  parallel  connected  single  cells  in  a  fuel  cell  stack 

cool  regarding  the  cooler  in  the  DMFC  system 

IE  regarding  anion  exchange  or  anion  exchanger 

in  physical  value  at  inlet 

FC  fuel  cell 

feed  fed  to  the  anode  compartment 

Me  methanol 

need  needed  for  stoichiometric  yield 

out  physical  value  at  outlet 

prod  produced  amount  of  the  physical  value 

reaction  lost  due  to  the  electrochemical  reaction  at  anode 

ret  implying  anion  retention 

total  total  physical  value 

xover  crossover  from  anode  to  cathode  side 
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